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^ . Abstract 

\ We study the structures in cross sections and R ratios concerning e + e~ annihilation 

reactions, for open-strangeness, open-charm and open-beauty hadron-pair production, from 
the perspective of a recently developed production formalism. Special attention is paid to 
non-resonant contributions. Known, tentative, and yet unknown resonances are spotted, 

\ with estimates for their masses and widths. 

> 

M 1 Introduction 



Several of the new enhancements in production cross sections share one common property, namely, 
they occur at — or just above — an important threshold. The most recent example, the J/ipcf) 
enhancement observed and baptized F(4140) by the CDF Collaboration [1], appears right above 
the J/ip4> threshold. The enhancement in the e + e~ — > A+A~ cross section reported by the 
Belle Collaboration [2] occurs right above the A+A~ threshold. The Y (4260) enhancement in the 
e + e~ — > J/vljn + 7T~ cross section, observed by the BaBar Collaboration [3], is right on top of the 
D* S D* S threshold. The X(3872) enhancement [4] in — > J/^/i ± 7r + 7r~ decay lies just above the 
DD* threshold. The new enhancement found by the BES Collaboration at 3.763 GeV [5] in the 
cross section of hadronic electron-positron annihilation reactions shows up just above the DD 
threshold. Therefore, we expect a common explanation for these phenomena. 

In the present work, we shall apply a recently developed formalism for the study of hadronic 
electron-positron annihilation reactions to determine common properties of the signal for open- 
strangeness, open-charm, and open-beauty hadron-pair production near threshold. The formal- 
ism is outlined in Sec. 2 and applied in Sec. 3. Conclusions are presented in Sec. 4. 
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2 Hadron-pair production in e e + annihilation 

In Ref. [6], we derived a relation between the production amplitude a^a — > i)s for e + e~ to two- 
meson annihilation reactions and the matrix elements of the meson-meson scattering amplitude 
T, resulting for the £-th partial wave in an expression of the form Tg(v — > i) 

a e (a -> i) oc g ai j l (pjr ) + ~ ^ (p v r Q ) g a uT e (i -> v) , (1) 

V 

where the g a i stand for the relative couplings of each of the two-meson systems i to a qq state 
of flavor a, and ji and hf~ are the spherical Bessel and the spherical Hankel function of the first 
kind, respectively. The interaction radius r$ represents the average distance at which the meson 
pair emerges from the interaction region. Furthermore, pi is the relative linear momentum in 
two-meson channel i. The matrix element T^iv — > i) of the scattering amplitude T describes 
transitions between channels v and i in inelastic meson-meson scattering. 

In the form presented in Eq. (1), the production amplitude consists of a non-resonant term, 
given by a Bessel function, and a sum of terms proportional to matrix elements of the scattering 
amplitude. Resonances reside in the latter matrix elements. Consequently, we have a clear sepa- 
ration of non-resonant and resonant terms for the production amplitude. We shall demonstrate 
in the following that such a separation is useful. 

Expression (1), which was derived in the framework of the Resonance-Spectrum Expansion 
(RSE) [7], is a direct link between quark dynamics and quantities observed in experiment. How- 
ever, in its present form the RSE scattering amplitude does not yet fully correspond to the 
physical reality, since it does not account well for the experimental observation [8] that two- 
meson channels only couple significantly to the quark- ant iquark propagator in a very limited 
range of energies close to threshold. From harmonic-oscillator (HO) confinement, which is the 
basis of the RSE, it is reasonably well understood why a specific channel does not survive at 
higher invariant masses. Namely, it is straightforward to derive [9, 10] that the number of dif- 
ferent channels to which the qq system couples grows rapidly for higher excitations. So for a 
given probability of quark-pair creation, the fraction that is left for a specific channel decreases 
accordingly. 

All terms of the production amplitude of Eq. (1) are proportional to three- meson couplings, 
given by g au . However, the three-meson couplings also enter the scattering amplitude T itself [6], 
that is, their squares. So at higher energies the terms containing the resonances fall off faster 
than the non-resonant term, leading to qq resonance peaks that fade away against background 
fluctuations of the non-resonant contribution in production cross sections. In Refs. [11, 12], we 
showed that signs of the ip{5S) and ip(AD) charmonium resonances can be observed just above 
threshold in the e + e~ — > A+A~ cross section measured by the Belle Collaboration [2], with 
tentative ip(6S) and ip(5D) signals, at about 400 MeV higher invariant masses, almost drowned 
in the background. 

By the use of Eq. (1) and the three-meson couplings of Refs. [9, 10], we may thus understand 
why it is hard to observe qq resonance peaks in scattering and especially production. However, 
from Ref. [8] we conclude that such a suppression may even be considerably larger, ranging from 
a factor 49 ± 25, determined via the process a (DD*) /a (D*D*), to a factor 124 ± 100, obtained 
from a(DD) ja (D*D*), thereby already accounting for the combinatorial factors PP : PV : 
VV = 1:4:7 and scaling invariant mass with the RSE HO frequency uj. Consequently, we must 
take one step back in order to figure out what extra ingredient is needed to comply with the 



2 



experimental evidence. Here, we focus on two-meson production reactions in e + e~ annihilation, 
which we assume to take place via the coupling of the qq propagator to the photon. 

In Refs. [9,10], a relation between intermesonic distances and the strength of the three-meson 
coupling for configuration space was derived, assuming HO confinement and 3 Pq quark-pair 
creation. However, since for HO distributions Fourier transforms are obtained by the simple 
substitution of coordinate x by linear momentum k, one readily gets the corresponding relation 
between the three-meson coupling and relative linear momentum. In Fig. 1, we depict the shape 
of the three-meson vertex intensity [10] for the cases to be considered here. In the same figure, 
we show that this shape, being an approximation itself, can be further approximated quite well 
by a Gaussian shape [13]. 




Figure 1: The shape of the three-meson vertex strength {MM \ V(k) \ M) as a function of the relative 
meson-meson linear momentum k (^i^), for V — > PP, V — > PV and V — ► VV (V = vector, P = pseu- 
doscalar) transitions [10], assuming HO confinement and 3 Po quark-pair creation, and its approximation 
by a Gaussian shape (^n^). The matrix element (MM \ V(k) \ M) and the linear-momentum variable k, 
are dimensionless here. 

For e~e + annihilation into hadrons in the charmonium region of invariant masses, we assume 
that such processes take place through the cc propagator, and, moreover, that the amplitudes 
are dominated by loops of open-charm mesons. The latter assumption implies that, even if the 
final state contains charmless hadrons, the line shape of its production cross section is to a large 
extent proportional to the line shape of open-charm production. 

The oscillator states must be decomposed into two- meson plane waves [14], which leaves us 
with the relative linear momentum p as the relevant parameter, rather than the total invariant 
mass. In the following we shall demonstrate that, for simple one-channel reactions, in the absence 
of inelasticity, like e.g. DD production in e + e~ annihilation far below the DD* threshold, we 
may reduce expression (1) for the resulting partial P-wave cross section o> to something of the 
form 

o P = 2A 2 4vra 2 r 2 

where, for the case of a reasonably narrow nearby resonance, such as e.g. if) (ID), we have sim- 
plified the part of Eq. (1) containing the scattering amplitude to a phase factor multiplying the 
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Breit-Wigner (BW) expression 



A(BW) 



r/2 

y/s-M + iT/2 



(3) 



here, \/s and M — iT/2 are the total invariant mass and the pole position of the resonance 
in the complex E plane, respectively. The dimensionless parameter A represents an overall 
normalization, while b adjusts the relative intensities of the resonant and the non-resonant term. 

In order to understand well the meaning of the phase in Eq. (2), we depict in Fig. 2 an 
elastic-scattering equivalent of Eq. (2). Of course, there is no elastic-unitarity requirement in 
production. Hence, unambiguous definitions of production phases cannot be made. Moreover, 
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Figure 2: Schematic Argand diagram for the phase motion of a two-meson production amplitude, for 
the case that no resonance is present (a), and for the case that there also is a resonance in the same 
invariant-mass interval (b). At threshold we assume zero phase. The arrows indicate the direction of 
increasing invariant mass. The angle <p, representing the angle of the phase motion of the non-resonant 
part of the amplitude at the central resonance mass, is defined in Eq. (2). 

disentangling non-resonant signals and resonance amplitudes is model dependent, even in the 
most naive approaches. Here, we propose to use our model [6] for which the scattering amplitude 
T contains the quark-antiquark propagator spectrum from the start, and for which some of the 
resulting resonances can be approximated by BW amplitudes. 

In Fig. 2a, we sketch the phase motion of a non-resonant signal, as an illustration. At 
threshold we assume a vanishing phase. The phase grows with increasing invariant mass, up to 
a certain maximum, after which it decreases. The maximum may be well over tp = 90°, but 
this does not necessarily imply that we are dealing with a resonance (cf. Levinson's theorem). 
The phase motion of a resonance, riding on top of the falling slope of the non-resonant signal, is 
indicated in Fig. 2b. While the phase of the non-resonant contribution moves on, the phase of the 
resonance makes a full 180° turn, to end up back on the slope of the non-resonant contribution 
at some higher invariant mass. The average phase ip of the non-resonant contribution in Eq. (2) 
is assumed to be a constant. 
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3 The parameter r 



The parameter r , which was introduced in Ref. [15] in order to allow for a quantity representing 
the average distance of quark-pair creation, requires renewed attention here. 

In Fig. 1 we have depicted the shape of the 3 Pq coupling as a function of the two-meson relative 
linear momentum p. We find that it vanishes at p — 0, rises to a maximum, and then drops 
fast for increasing values of p. This behavior of the three-meson vertex intensity as a function 
of two-meson relative linear momentum agrees with the experimental observation obtained with 
the SLAC/LBL magnetic detector at SPEAR [8]. The results suggest that at higher invariant 
masses the branching fraction of a specific two-meson channel vanishes rapidly. However, the 
rate at which the branching fractions decrease still remains to be studied. 

As an example, let us concentrate on the reaction e + e~ — > BB. In the RSE, the BB threshold 
at 10.56 GeV comes out far above the HO ground state at 9.72 GeV and the first radial excitation 
at 10.10 GeV, and also above the second radial excitation at 10.48 GeV. Note that the radial 
level spacings are equal to 2u = 0.38 GeV in the RSE. 

The three-meson vertex intensities for these HO levels can be determined with the formalism 
developed in Ref. [10]. The total BB branching fraction is found to be roughly 2.8% at the HO 
ground state, but only about 0.4% at the second radial HO excitation. The main competition for 
BB formation above the BB threshold stems from the BB* channel. This competing channel 
will start to dominate for invariant masses closer to the BB* threshold at 10.60 GeV. Above the 
BB* threshold, we suppose that BB formation rapidly vanishes. 

The above implies that at the BB* threshold the coupling of BB to the bb propagator must 
have decreased significantly. This is indeed confirmed by Fig. 3, in which we show data for the 
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Figure 3: Experimental data for the process e + e~ — > bb obtained by the BaBar Collaboration in 
Ref. [16]. The vertical lines indicate the BB* and B*B* thresholds, as indicated in the figure. The 
eye-guiding lines reflect our interpretation of the data, and do not represent fits to the data. 



process e + e~ — > bb, measured and analysed by the BaBar Collaboration [16]. As also remarked 
in their paper, the large statistics and the small energy steps of the scan make it possible to 
clearly observe the two dips at the opening of the thresholds corresponding to the BB* + BB* 
and B*B* channels. 

We have cut off the huge peak at 10.58 GeV, in order to concentrate better on the details of 
the other two enhancements, at 10.63 GeV and 10.69 GeV, respectively. Near the BB* threshold, 
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we thus observe that the BB signal rapidly vanishes for increasing invariant mass, whereas the 
BB* signal behaves the opposite way, i.e., it grows fast for increasing invariant mass just above 
threshold. At the B*B* threshold, this phenomenon is repeated, now with respect to the BB* 
signal. 

So we conclude that BB production in e + e~ annihilation can mainly be observed within the 
invariant-mass window delimited by the BB and BB* thresholds. Similarly, BB* production 
has an equally wide window formed by the BB* and B*B* thresholds. The window for B*B* 
production is somewhat wider, since the next threshold concerns the B S B S channel, which lies 
considerably higher. Therefore, the enhancement peaking at about 10.69 GeV is broader than 
the ones at 10.58 and 10.63 GeV. 

From Fig. 1, we see that the bb —> BB form factor closes its window at about k = 2.5. As 
a result of the above discussion and considering the structure of the production amplitude in 
Eq. (1), this implies that the dimensionless quantity defined by 

p 2 BB (at M (BB) = Mb* + M B ) r 2 0BB = 

= {I W~ + ~ Ml) rl BB = \ (M B , - M 2 ) {l + j^^} r 2 , B§ (4) 

must be equal to about 2.5 2 ~ (0.5 GeV im/hc) 2 . 

Owing to quadratic meson-mass relations, the quantity My — Mp (V = vector, P = pseu- 
doscalar meson) seems to be almost a constant of nature for the cases we shall consider here, 
namely 0.488 ± 0.007 for V = B* and P = B, 0.550 ± 0.001 for neutral and 0.546 ± 0.002 for 
charged mesons, and 0.5513 ± 0.0005 for K^*^. Furthermore, for heavy quarks, the expres- 
sion between braces on the righthand side of Eq. (4) approximately equals 2. This suggests that 
we could set out by assuming that the whole quantity in Eq. (4) is a constant. In the following we 
shall show that indeed the constant p P p (at M = My + M P ) r 0PP = 2.5 leads to good results 
for P = D,B (V = D*, B*), and even for P = A c {V = S c ), though for P = K (V = K*) we 
need a somewhat smaller value for this constant. 

In view of the above, we find for the tq parameter in the cc and bb sectors considerably larger 
values than we had assumed in earlier work. Although this is not entirely unwelcome, as it may 
help cure the somewhat too small widths of cc and bb resonances we obtained before, it seems 
to break, to some extent, the rigorous flavor independence we advertised in the past. A more 
detailed discussion of how to precisely interpret ro is required. The problems lies in the fact 
that the valence particles involved in the decay process change their character and in general 
their masses, that is, an original qq pair transforms, in the course of the decay process, into 
a pair of mesons. In a coordinate-space description, which was the original formulation of the 
model [15, 17], it turned out necessary to assume that the r coordinate of the valence qq pair is 
the same as for the emerging meson pair, in order to avoid nonlocalities in the equations. This 
approximation is reasonable when both the decaying meson and the decay products are composed 
of light quarks, so that all mesons have comparable sizes. However, when the original meson is 
a heavy quarkonium, which is a much more compact system especially in the bottomonium 
case, the main decay products are mesons with open beauty (or charm), which are clearly much 
more extended objects owing to the presence of a light quark. In these situations, the local 
approximation is obviously not ideal, but the tendency for ro will be to adjust itself chiefly to the 
kinematics of the outgoing meson pair, especially when data over a wide p range are analysed. 
Possible improvements will be left to future work. One alternative is the employment of more 
extended transition potentials [10], so as to allow for more detail in the r dependence of the 
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wave functions. Another possibility is to work in a covariant momentum-space framework [18], 
in order to be able to deal with nonlocalities. 

In view of these considerations, it seems most appropriate to consider ro, which appears in the 
Bessel and Hankel functions of Eq. (1), a parameter exclusively governed by the meson-meson 
dynamics and kinematics, at least for the purpose of the present analysis. The mild deviations 
from exact flavor independence of the meson spectra, resulting from this approach, are irrelevant 
for this work, and would generally have only minor effects on the central resonance positions 
predicted by the RSE. 

In the following we compare e + e~ annihilation data above the open-strangeness threshold for 
K + K~ production [19,20], just above the open-charm threshold for inclusive hadron production 
[5], just above the open-beauty threshold for inclusive hadron production [21] and exclusive bb 
production [16], and just above the open-charm baryon-pair threshold for A+A~ production [2]. 

The production cross sections of hadron pairs are expressed in nanobarns in the above exper- 
imental analyses. Hence, we must scale expression (2) with the size of the available window for 
production and by the appropriate conversion factors. Thereto, by the use of the quantity given 
in Eq. (4) and its equivalent for the other processes studied here, we perform in expression (2) 
the substitution 

i ra x 10 7 nb/fm 2 

AVg^A' 2 , (5) 

y2.5 2 + (m h r f - m h r 

in order to facilitate direct comparison of Eq. (2) with experimental results. For m/, given by 
?77, £>, rriB, ttia c , arid m,K, we next adjust the shape in Eq. (2) to the shapes of the experimental 
cross sections, by choosing for A' the values 0.95, 1.03, 1.02, and 0.49 (only K + K~ has been 
measured), respectively. 



3.1 The reaction e + e — ► DD 

The reaction of e + e~ annihilation into open-charm pairs can be observed at and above the DD 
threshold. We assume here that the reaction takes place via a photon and the cc propagator, 
through the creation of a light qq pair. However, as we have discussed previously, many competing 
configurations may be formed, increasing in number for higher invariant masses. Furthermore, it 
seems we may conclude from experiment that stable open-charm hadrons have more probability 
to be formed near threshold, where kinetic energy is almost zero. Hence, if it were not for phase 
space and the centrifugal barrier, DD pairs would be produced most likely just above threshold. 

For total invariant mass below, but close to the DD* threshold, we assume that the proba- 
bility of DD creation is already reduced because of the non- vanishing probability of creating a 
virtual DD* pair. Just above the DD* threshold, DD creation decreases rapidly. The expected 
corresponding non-resonant contribution of the production amplitude should thus exhibit this 
feature. In this respect, an important observation was published by the BES Collaboration in 
Ref. [5]. To our knowledge, the BES Collaboration was the first to discover that the ^(3770) cross 
section is built up by two different amplitudes, viz. a relatively broad signal and a rather nar- 
row cc resonance. For the narrow resonance, which probably corresponds to the well-established 
if) (ID) (3770), the BES Collaboration measured a central resonance position of 3781.0 ± 1.3 ± 0.5 
MeV and a width of 19.3 ± 3.1 ± 0.1 MeV (their solution 2). If the latter parameters are indeed 
confirmed, it would be yet another observation of a quark-antiquark resonance width that is very 
different from the world average (83.9 ± 2.4 MeV [22] in this case), after a similar result was 
obtained by the BaBar Collaboration in Ref. [16], for bb resonances. Concerning the broader 
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structure, the BES Collaboration indicated, for their solution 2, a central resonance position of 
3762.6 ± 11.8 ± 0.5 MeV and a width of 49.9 ± 32.1 ± 0.1 MeV. The signal significance for the 
new enhancement is 7.6a (solution 2). It is has been explained as a possible diresonance [23] or 
heavy molecular state [24]. 

Furthermore, in the latter BES publication, the existence of conflicting results with respect 
to the branching fraction for non-DD hadronic decays of the ip(\D)(3770) was emphasized. On 
the one hand, the total branching fraction for exclusive non-DD modes has been measured to be 
less than 2% [25,26]. But on the other hand, for inclusive non-DD decay modes, values of about 
15% have been found [27,28]. According to the BES Collaboration, this apparent discrepancy 
may partially be due to the assumption that the line shape above the DD threshold is the result 
of one simple resonance. Here, we shall show that the broader structure is most likely caused by 
the non-resonant contribution to the production amplitude given in Eq. (2), thus lending further 
support to the idea that the ijj(lD)(3770) enhancement consists of two distinct signals. 

We mentioned in the above the recent findings of the BaBar Collaboration with respect to 
the resonance parameters of bb states. The corresponding measured line shapes actually show a 
similar behavior, namely a broad structure with a narrow resonance on top. Nevertheless, the 
BaBar Collaboration did not conclude that the broad structures are due to additional resonances 
in the bb sector. As we shall demonstrate in the following, non-resonant contributions can indeed 
be large, sometimes even much larger than the qq resonances in their vicinity. 

If we assume the above-mentioned value of 2.5 for the quantity in Eq. (4), we get r = 0.962 fm. 
A comparison of the resulting theoretical curve, by the use of Eq. (2), together with the data 
of the BES Collaboration, is shown in Fig. 4. The BW parameters of the ^(3770) used for the 




Figure 4: Comparison of experimental data obtained by the BES Collaboration [5] and the cross 
section resulting from Eq. (2) (^»^) for tq = 0.962 fm and b = 1.33. The non-resonant contribution 
^) dominates for larger relative DD linear momentum. The BW parameters of the cc resonance are 
M(ip(lD)) = 3.781 GeV and r(^(LD)) = 17 MeV. 

theoretical curve of Fig. 4, are collected in Table. 1 and compared to the corresponding PDG 
parameters of those of the BES Collaboration. 
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source 


Mass (GeV) 


Width (MeV) 


phase (rad) 


Fig. 4 


3.781 


17 


0.75tt 


BES [5] 


3.7810 ± 1.3 ±0.5 


19.3 ±3.1 ±0.1 (Oi 


>8 ± 1.86 ± 0.03)tt 


PDG [22] 


3.7700 ±0.0024 


83.9 ±2.4 





Table 1: Parameters used for the solid line in Fig. 4, compared to the corresponding parameters obtained 
by the BES Collaboration [5], and to world averages [22]. 

3.2 The reaction e + e~ BB 

The higher excitations of the beautonium vector states, discovered more than two decades ago, 
are still today a puzzling topic of intensive research. In Refs. [29] and [30], the CUSB and 
CLEO Collaborations, respectively, displayed the first results for the invariant-mass spectra of 
the R (ohad/o^j) ra ti° above the open-beauty threshold. 

The data of Ref. [29] were observed with the CUSB calorimetric detector operating at CESR 
(Cornell). The experimental analysis resulted in evidence for structures at 10577.4 ± 1 MeV, 
10845 ± 20 MeV, and 11.02 ± 0.03 GeV, with total hadronic widths of 25 ± 2.5 MeV, 110 ± 15 
MeV, and 90 ±20 MeV, respectively. Structures at about 10.68 and 11.2 GeV were not included 
in the analysis of the CUSB Collaboration. 

The data of Ref. [30] were obtained from the CLEO magnetic detector, also operating at 
CESR. The experimental analysis resulted in evidence for structures at 10577.5 ± 0.7 ± 4 MeV, 
10684 ± 10 ± 8 MeV, 10868 ± 6 ± 5 MeV, and 11019 ± 5 ± 5 MeV, with total hadronic widths of 
20 ± 2 ± 4 MeV, 131 ± 27 ± 23 MeV, 112 ± 17 ± 23 MeV, and 61 ± 13 ± 22 MeV, respectively. A 
structure at about 11.2 GeV was not included in the analysis of the CLEO Collaboration. 

The enhancement at 10.580 GeV was extensively studied in a more recent publication [21] of 
the BaBar Collaboration. Data were collected with the BaBar detector at the PEP-II storage 
ring of Stanford Linear Accelerator Center. The experimental analysis yielded 10579.3 ± 0.4 ± 1.2 
MeV and 20.7 ± 1.6 ± 2.5 MeV for the central mass and the total (hadronic) width, respectively, 
which in in fair agreement with the above results of the CUSB and CLEO Collaborations. 

A very recent release of the BaBar Collaboration [16] is devoted to two more resonances in 
the bb spectrum, also using data obtained by the BaBar detector at the PEP-II storage ring, 
resulting in the BW parameters 10876 ± 2 MeV (mass) and 43 ± 4 (width) for the T(10860), 
and 10996 ± 2 MeV (mass) and 37 ± 3 (width) for the T(11020), which differ substantially, in 
particular for the widths, from the above results of the CUSB and CLEO Collaborations, and 
also from the world-average values [22] 10865 ± 8 MeV (mass) and 110 ± 13 (width) for the 
T(10860) and 11019 ± 8 MeV (mass) and 79 ± 16 (width) for the T(11020). Of course, such 
discrepancies call for further study. However, in this work we shall only focus on the reaction 
e + e~ — > BB for total invariant masses below the BB* threshold. 

The data of the BaBar Collaboration of Ref. [21] consist of three different scans of the visible 
cross section for bb production, with a background of other hadronic processes in e + e~ annihila- 
tion. The three scans differ in the amount of probable background, while the first scan also has 
its peak shifted. Here, we subtract backgrounds of 0.82 nb, 0.91 nb, and 0.99 nb, respectively, 
from the data of scans I, II and III. Moreover, we shift the data of scan I by 4 MeV. The shape 
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of the resulting signal of 17 data points above the BB threshold is compared to the data of the 
BaBar Collaboration [16]. For the latter data, which show the Rb ratio for all e + e~ annihilation 
processes containing b quarks, we assume a background of Rb = 0.1, to account for those reac- 
tions that do not contain BB pairs. Fig. 5a shows the resulting line shape of the enhancement 
at 10.58 GeV. 




10.56 10.58 10.60 10.56 10.58 10.60 

M(GeV) M(GeV) 



Figure 5: Comparison of experimental data taken from (a) Refs. [21] (•) and [16] (v)> an d (b) Refs. [29] 
(y) and [30] (■), with the line shape resulting from Eq. (2) for tq = 1.014 fm. The vertical line in (b) 
indicates the position of the BB* threshold. 

The data of the CUSB [29] and CLEO [30] Collaborations are shown in Fig. 5b. The former 
data are as measured, but subtracted by invisible ) = 2.3 in order to account for background, 
as suggested by CUSB themselves. Similarly, the latter data are as measured, but subtracted 
by R = 4.57 to account for background, as suggested by CLEO. We have not put the CUSB 
and CLEO data in one figure together with the BaBar ones, since it is not clear to us how to 
normalize each separate data set so as to make sense out of a joint presentation. Moreover, in 
Fig. 5b the CUSB and CLEO data are treated as two different sets, each compared separately 
to the shape resulting from Eq. (2). 

The theoretical curves in Fig. 5 correspond to the non- resonant contribution in Eq. (2). In 
Ref. [22], this enhancement is classified as a bb resonance, under T(10580). However, in view of 
our results, we do not believe this enhancement to represent a resonance. In fact, we expect [15] 
the corresponding 4 3 S*i bb resonance to lie at roughly 10.77 GeV. A better candidate for this 
bb resonance is the state found by the CLEO Collaboration [30] at 10684 ± 10 ± 8 MeV. It is 
unclear to us why this resonance, observed back in 1985, was classified as a bbg hybrid state by 
the CLEO Collaboration, and, subsequently, not even mentioned in the PDG tables. 

The enhancement at 10.58 GeV suggests an accumulation of BB pairs in this invariant- mass 
region. Therefore, a description in terms of a wave function with a dominant BB component 
appears to be more adequate than assuming a pole in the scattering amplitude due to a supposed 
underlying bb state. 

3.3 The reaction e + e~ —*■ A+A~ 

In Ref. [2], the Belle Collaboration announced the observation of a near-threshold enhancement, 
by studying the e + e~ — > A+A~ cross section. The experimental analysis resulted in a mass and 
width for this enhancement of M = (4634lf)(stat.)i|(sys.) MeV and T tot = 92^(stat.)±2?(sys.) 
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MeV, respectively [2], with a significance of 8.8 a. An intriguing aspect of this experimental 
observation is that the main signal lies close to the A+A~ threshold, making an understanding 
of this structure a highly topical issue. 

An initial study of the Belle cross section was done in Ref. [11]. Here, we apply Eq. (2), but 
assume an S rather than a P wave, by choosing the appropriate Bessel and Hankel functions. The 
result is shown in Fig. 6, where, besides reproducing the enhancement at threshold, we also spot 




4.&0 ' 4?T0 4.&0 4.90 

M(A+A-)(GeV) 

Figure 6: Comparison of the experimental cross section for the reaction e + e~ — > A+A~, obtained by 
the Belle Collaboration [2], and the line shape resulting from the S- wave equivalent of Eq. (2) 
for ro = 0.962 fm (see also Ref. [13]). The non-resonant contribution ( «»■,) is relatively large. The 
BW parameters for the two cc resonances [11] are M (ip(5S)) = 4.784 GeV, T (ip(5S)) = 55 MeV, and 
M (ip{4D)) = 4.871 GeV, T (V>(4L>)) = 60 MeV. 

two new charmonium states, identified as the ip{5S) and ip(4D). We get 4.784 GeV and 4.871 
GeV for the central resonance masses, and 55 MeV and 60 MeV for the widths, respectively. 
Furthermore, the curves in Fig. 6 clearly show that the presence of the ip{5S) and ip(4D) cc 
resonances strongly affects the production of A+ A~ near threshold. The latter shape will be 
further influenced by the ip(3D) cc resonance around 4.53 GeV [11,15,31], which still needs 
confirmation. However, this state cannot be parametrized with a simple BW structure, which 
does not have the right threshold behavior. Nevertheless, the full RSE amplitude, without the 
approximation of resonances by BWs as employed in the present analysis, may be used for the 
study of such effects. This will be done in future work. 

Modeling the the A+A~ enhancement could of course be done by considering a full wave 
function with all possible components, viz. cc states, charmed-meson pairs, and charmed-baryon 
pairs. Such a wave function will have a large A+A~ component under the enhancement. Never- 
theless, it will not give rise to a resonance pole in the full coupled- channel scattering amplitude, 
unless, accidentally, there is a dynamically generated resonance pole in this invariant-mass re- 
gion. Now, in view of the large non-resonant contribution to the A+A~ enhancement, we do not 
consider the occurrence here of a dynamically generated resonance very likely, which is in line 
with a similar conclusion of D. V. Bugg [13]. In our opinion, it certainly does not represent a 
state of the J PC = 1 cc spectrum as advocated in Refs. [32,33]. 
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3.4 The reaction e + e — ► K + K 

Cross sections measured for the reaction e + e~ — > K + K~ date back to the 1970s. Here, we shall 
compare Eq. (2) to data from the ORSAY-DCI-DM1 experiment for invariant masses ranging 
from 1.4 GeV to 2.06 GeV, published in 1980 [20], data measured by the NOVOSIBIRSK-OLYA 
experiment in the energy range 1.017-1.4 GeV, published in 1981 [19], more recent data from the 
CMD-2 Collaboration determining the 0(1020) mass and width [34], published in 1995, and very 
recently published (2007) data obtained with the SND detector in the VEPP-2M experiment, for 
invariant masses ranging from 1.04 GeV to 1.38 GeV [35]. The results of Eq. (2) for tq = 0.56 
fm are shown in Fig. 7. The BW parameters of the four structures of sufficient significance in 
our analysis of the data are collected in Table 2. 




2.0 
1.5 
1.0 
0.5 



0.2 0.4 0.6 

p(K+K-)(GeV) 



? 




(c) 


CO 

O 


M 0(1020) 




b 













0.12 0.14 0.16 

p(K + K-)(GeV) 



-1 



1 

^b 

£ i 

o 

tH 






(b) 




/ 


(t 





Xi 
b 

10 



0.2 0.4 0.6 0.8 

p(K+K-)(GeV) 

KK* K*K* 







(d) 




X(1250) 
















X(1590) 








X(1830) 















1 i 

0.4 


0.6 


0.8 



p{K + K-){GeV) 



Figure 7: Comparison of combined experimental data taken from Refs. [20] (•), [19] (•), [34] (■), 
and [35] (•), with the line shape resulting from Eq. (2) (^ » ^), using the BW parameters of Table 2. 
On the horizontal axis we mark out the measured relative kaon momentum in the K + K~ CM system, 
and on the vertical axis the corresponding cross section, (a): Non- Resonant contribution (^n^) and 
sum of non-resonant and p meson contributions, (b): Including also the 0(1020) contribution, (c): The 
invariant-mass region of 0(1020). (d): Contributions of further resonances. 



Although we expected the KK* channel to be the main competitor for the production of 
charged-kaon pairs, we do not even observe any relevant decrease in the K + K~ cross section at 
the KK* threshold. On the contrary, the non-resonant contribution to the K + K~ line shape 
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is maximum at the KK* threshold. Indeed, from a detailed study performed by the BaBar 
Collaboration in Ref. [36] we learn that, while the individual isoscalar and isovector cross sections 
for e + e~ — > K*K are comparable in magnitude to those for e + e~ —* K + K~, the total cross 
section for e + e~ — > K + K~tt° is very small, hinting at destructive interference between isoscalar 
and isovector contributions. Moreover, several other channels compete with the production of 
K + K~ , like rju, rjcj), rj'u, rj'cj), irp in the isoscalar sector, and rjp, rj'p, pp in the isovector sector. 
Several of these channels are in reality smeared out, as they involve the broad p resonance. 
Consequently, their effect is spread out over a relatively large interval of invariant K + K~ masses. 
The strongest and sharpest channel is probably 7/0, with threshold at almost 2.0 GeV. Hence, 
it is not a real surprise that our simple rule of thumb for the quantity defined in Eq. (4) does 
not work in this case. For the non-resonant contribution to the theoretical cross section shown 
in Fig. 7, we have chosen tq = 0.56 fm, which is equal to the kaon charge radius. 

A further consequence of the apparent absence of channels competing with K + K~ production 
is that it allows us to analyse the experimental cross section well above threshold by using Eq. (2), 
thereby assuming that no relevant inelasticity occurs. 

In Fig. 7a we show the two dominant "background" contributions, viz. the non-resonant 
term of Eq. (2) and the contribution of the p resonance. The 0(1020) resonance is included 
in Fig. 7b, using the PDG resonance parameters [22]. In Fig. 7c only the 0(1020) signal is 
shown. The theoretical cross section is adjusted to the data in Fig. 7d by introducing three 
resonant structures, namely (also see Table 2) X(1250) (M = 1.25 GeV, T = 60 MeV), X(1590) 
(M = 1.59 GeV, T = 60 MeV), and X(1830) (M = 1.83 GeV, T = 80 MeV). Additionally, we 
display in Fig. 7d the contribution of the non-resonant term of Eq. (2) 



piek 


Mass (GeV) 


Width (MeV) 


0(1020) 


1.01946 


4.26 


X(1250) 


1.25 


60 


X(1590) 


1.59 


60 


X(1800) 


1.83 


80 



Table 2: BW parameters used for a fit to the data shown in Fig. 7. The 0(1020) parameters are taken 
from Ref. [22]. 

One observes from Fig. 7 that data for the reaction e + e~ — > K + K~ roughly follow the 
tail of the non-resonant term of Eq. (2) for K + K~ momenta larger than about 0.6 GeV, which 
corresponds to a total invariant mass of 1.56 GeV. Below this mass, the 0(1020) and p resonances 
dominate the cross section. The 0(1020) line shape is well reproduced by our theoretical curve, 
which uses the PDG [22] values for its mass and width. 

The recent data of Ref. [35] are in conflict with the almost three decades older data of 
Ref. [19]. Our theoretical curve agrees better with the former data for p(K + K~) « 0.19 GeV, 
which corresponds to an invariant mass of 1.07 GeV. In particular, three subsequent data points 
of Ref. [19] do not fit in the amplitude of Eq. (2). However, in the invariant-mass region of 1.17— 
1.35 GeV (p(K + K~) ~ 0.30-0.45 GeV), we notice that the data of Ref. [35] are substantially 
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larger than those of Ref. [19] and, moreover, that some of the structure has flattened out. We 
should have no difficulty in fitting the data of Ref. [35] , as it would only change the parameters 
of the X(1250) resonance. Here, we have opted for reproducing the data of Ref. [19] in the 
invariant-mass region of 1.17-1.35 GeV. 

4 Conclusions 

We have shown that our recently developed amplitude for hadron-pair production in electron- 
positron annihilation describes quite well the corresponding experimental data for open strange- 
ness, open charm, and open beauty. Moreover, the formalism allows for a well-defined separation 
in resonant and non-resonant contributions. By choosing the value 2.5 for the quantity that 
defines the effectively available window for production, given in Eq. (4), we tune the only free 
parameter r controlling the non-resonant contributions to data for DD, BB, A+A~, and K + K~ . 
As a result, we obtain for r a value close to 1 fm for the former three process, whereas for the 
latter one the kaon charge radius turns out to work well. Finally, it follows from our approach that 
the size of the window rather accurately defines the total amount of non-resonant hadron-pair 
production. 

In the DD case, we find, besides a substantial non-resonant contribution, a narrow ^(3770) 
resonance, with a central mass of 3.781 GeV and a width of 17 MeV, both well within the 
error bars suggested by the BES Collaboration. As for the bottomonium sector, the T( 10580) 
enhancement in BB appears to contain exclusively non-resonant contributions, so that it should 
not be associated with a bb state. A wave function with a large BB component is likely to be a 
more appropriate description of this enhancement. 

For the cross section of A+A~ production, one needs to assume, besides the non-resonant 
enhancement near threshold, two new cc resonances, which are naturally interpreted as the 
ip(5S) and the ip(4D). We obtain 4.784 GeV and 4.871 GeV for their central resonance masses, 
and 55 MeV and 60 MeV for their widths, respectively. The non-resonant contribution to K + K~ 
production is rather modest. Apart from a very accurate description of the 0(1020), we find three 
further enhancements that can be fitted as BWs, namely the X(1250), X(1590), and AT(1830), 
with their central masses at 1.25 GeV, 1.59 GeV, and 1.83 GeV, respectively, and widths of 60-80 
MeV. 
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